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Abstract 

A majority of the known Colias species (Lepidoptera: Pieridae, Coliadinae) occur in the mountainous 
regions of Centtal-Asia, vast ateas that are hard to access, rendering the knowledge of many species limited 
due to the lack of extensive sampling. Two gene tegions, the mitochondrial COI 'barcode' region and the 
nuclear ribosomal protein RpS2 gene region wete used fot exploting the utility of these DNA markers 
fot species identification. A comprehensive sampling of COI barcodes fot Centtal Asian Colias buttetflies 
showed that the barcodes facilitated identification of most of the included species. Phylogenetic recon- 
struction based on patsimony and Neighbour-Joining recovered most species as monophyletic entities. 
Fot the RpS2 gene region species-specific sequences were registered for some of the included Colias spp. 
Nevertheless, this gene region was not deemed useful as additional molecular 'barcode'. A parsimony 
analysis of the combined COI and RpS2 data did not support the current subgenetic classification based 
on motphological chatactetistics. 
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Introduction 

The use of a standardized gene region, i.e. a 650 bp fragment of the 5'-region of the 
mitochondrial cytochrome c oxidase subunit I (hereafter COI), as a DNA barcode 
(Hebert et al. 2003), to facilitate identification of biological specimens, as well as for 
calling attention to possible new species has generated a steadily increasing number of 
DNA barcoding studies of invertebrates (Taylor and Harris 2012), and particularly of 
Lepidoptera (see www.lepbarcoding.org). While the utility of DNA barcoding as an 
investigative tool has gained much support, there still remain a number of problems 
related to the use of a single DNA sequence as a taxon barcode. Several studies on 
Lepidoptera have shown that species may be polymorphic and/or share haplotypes 
(Nice et al. 2002, Wahlberg et al. 2003, Elias et al. 2007, Schmidt and Sperling 2008), 
so that identifications may become less reliable. Additionally, it has been shown that 
incomplete lineage sorting or mitochondrial introgression could obscure the delimita- 
tion of closely related taxa (Tautz et al. 2003, Zakharov et al. 2009). Using one or a 
few specimens as representatives of a species indeed provides us with little information 
about their intraspecific variation, particularly for widely distributed species (e.g. Funk 
and Omland 2003, Seberg et al. 2003, Sperling 2003). 



The genus Colias 

The butterfly genus Colias Fabricius, 1807 is a genus of the family Pieridae (subfamily 
Coliadinae), comprising about 85 species. Most of its species have a limited distribution 
in the Arctic and Alpine regions of the Holarctic realm, but two species occur in the 
Afrotropical and seven are known from the Neotropical regions (Verhulst 2000). A few 
species are widely distributed and common, such as the Palaearctic C. erate (Esper, 1805) 
and C. croceus (Geoffroy, 1785), and the Nearctic C. eurytheme Boisduval, 1852 and C. 
philodice Godart, 1819. As a consequence, these taxa are frequently used in ethological, 
ecological and genetic research (e.g. Pollock et al. 1998, Wang and Porter 2004, Porter 
and Levin 2010). Colias erate and C. croceus are a species pair where only typical speci- 
mens can be reliably distinguished morphologically, and members of these species are 
known to frequently hybridize (e.g. Dinca et al. 20 1 1 and references therein) . Lukhtanov 
et al. (2009) indicated that mitochondrial introgression was a likely explanation for the 
shared barcodes they registered between these sympatric taxa. The Nearctic taxa C eury- 
theme and C. philodice are broadly sympatric sister species that hybridize frequently and 
that likely share a significant portion of their genomes through introgression (e.g. Wang 
and Porter 2004, Porter and Levin 2010). Verhulst (2000) illustrated hybrid individuals 
of six species of Colias from the Palaearctic region, including C. croceus. 

The Central Asian mountainous regions harbour nearly half of all Colias species. 
The distribution, ecology and taxonomy are still incompletely documented for most of 
these species, mainly due to their remote occurrences (Verhulst 2000). Central Asian 
Colias species occurring in remote mountainous areas that are hard to access have been 
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far less studied than their North American or European congeners. An important part 
of the older material that exists in museum collections worldwide (e.g. from Tibet) 
originates from early collecting expeditions in the late 19 th and early 20 th centuries. 
Important material was, however, also collected within the former Soviet Union dur- 
ing 20 th century. Fieldwork in Central-Asia has subsequently become less complicated, 
and thus new material is again available for research. As a result of this, new species 
such as Colias aegidii Verhulst, 1990 and Colias adelaidae Verhulst, 1991, have been 
described, as well as a number of new subspecies. Despite an increasing research effort 
on Central Asian Colias species there are as yet no published studies on their phyloge- 
netic relationships. 

The first contribution to the species classification of Colias was given by Berger 
(1986), who used a few morphological characters to establish a comprehensive subge- 
neric classification, comprising the subgenera Colias Fabricius, 1807, Neocolias Berger, 
1986, Eucolias Berger, 1986, Eriocolias Watson, 1895, Palaeocolias Berger, 1986, Si- 
milicolias Berger, 1986, Scalidoneura Butler, 1869 and Paracolias Berger, 1986. Later, 
Ferris (1993) used 84, mainly morphological, characters to reconstruct a phylogeny 
of all North American Colias species known at that time, which was the first spe- 
cies phylogeny within the genus Colias. The first contribution to the knowledge of 
the molecular phylogenetic relationships of the North American Colias species was 
made by Pollock et al. (1998), who studied a number of Colias species using a 333 
bp sequence fragment of the mtDNA COI gene. They found some small differences 
between species classified in the subgenera Neocolias and Eriocolias, thus supporting 
Berger's (1986) separation of Neocolias from Eriocolias. Pollock etal. (1998) also noted 
that even though Colias is a speciose genus, this was not mirrored in the COI sequence 
diversity. Wheat and Watt (2008) studied the molecular phylogenetic relationships of 
North American Colias taxa using mitochondrial gene sequences (ribosomal 12S and 
16S rRNA, Leu2 and Val tRNA and COI + II). Their results showed that the COI 
sequences only allowed identification of some of the taxa supported by the full data 
set used in their study. The results of their study further suggested that species radia- 
tions within Colias are comparatively young as compared with those of related pierid 
butterflies, since molecular divergences among species were small. Based on molecular 
data Brunton (1998) studied the phylogenetic relationships of the 12 Colias species 
occurring in Europe. He recovered three monophyletic groups largely corresponding 
to geographical distributions. He concluded that the Scandinavian species appeared to 
be the oldest in Europe, sharing a common ancestor with Colias species from the USA. 
According to Brunton (1998) the European Colias species radiated from Scandinavia 
to the rest of Europe forming an eastern clade and a western clade. As with Pollock 
et al. (1998), the results did not agree with Berger's (1986) subgeneric classification. 

The aim of the present study was to test the usefulness of COI barcodes for species 
identification of a broad representation of Central Asian Colias species, including nine 
Colias species overlapping with Lukhtanov et al.'s (2009) study, and 19 species not 
previously barcoded. In addition, we wanted to elucidate the informativeness of the 
RpS2 gene region that Wahlberg and Wheat (2008) found informative for lepidopter- 
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an phylogenetic relationships. We tested the nuclear ribosomal protein gene RpS2 as a 
potential complementary barcode region for Colias and for use in a combined analysis 
with COI for testing the current subgeneric classification of the species in the present 
study. We also contrasted our COI barcodes against a larger set of COI barcodes of 
Colias taxa available from GenBank (GB). 



Materials and methods 

Study area and taxon sampling 

This study includes material from the mountain regions of Kirgizistan, Tadzhikistan, 
northern Afghanistan, northern Pakistan and India (e.g. mountain ranges Tian Shan, 
Hindu Kush, Karakorum, Himalaya) and the mountain regions in the Chinese provinc- 
es Qinghai, Gansu, Sichuan, Yunnan and the autonomous regions Tibet and Xinjiang 
Uygur. The Colias fauna of these Central Asian regions comprises about 34 species (Ver- 
hulst 2000) while the species number for Central Asia in broad sense is over 40 species. 

The taxon sampling aimed to cover as many of the Colias species from this area as 
possible. Additionally, a few Colias species occurring in adjacent territories (e.g. Burya- 
tia) were also available for molecular study. Whenever possible, several individuals of 
each species were analysed to assess intraspecific variation. The available specimens 
used for molecular study consisted of a total of 56 adult specimens covering 27 species 
of Central Asian Colias and two Colias species from adjacent territories (Table 1). The 
specimens are preserved as DNA voucher specimens and labelled accordingly, to be 
deposited in the collections of the Zoological Museum of Finnish Museum of Natural 
History, Helsinki, Finland (MZH) (DNA voucher specimens MZHJL1-JL71). Spe- 
cies identifications were verified by JL based on easily recognizable diagnostic char- 
acters using the monograph by Verhulst (2000), while the taxonomy is according to 
Grieshuber and Lamas (2007). Additionally, we used 35 COI barcode sequences (17 
species) of Palaearctic Colias species obtained from GB, as listed in Table 2. 



Laboratory methods 

Total genomic DNA was extracted form 2-5 legs of dried, pinned butterfly specimens 
using NucleoSpin® Tissue Kit (Machery-Nagel), according to manufacturer's proto- 
cols, and resuspended in 50 ul ultrapure water. 

The primer pair LCO-1490 (5'-GGTCAACAAATCATAAAGATATTGG-3') 
and HCO-2198 (5'-TAAACTTCAGGGTGACCAAAAAATCA-3') (Folmer et al. 
1994) was used to amplify a ca. 650 bp fragment of the mitochondrial COI gene. The 
polymerase chain reactions (PCR) were done under the following parameters: initial 
heating 95 °C for 2 min, following 30 cycles of 94 °C for 30 s, 49 °C for 30 s and 72 
°C for 2 min, followed by a final extension of 72 °C for 7 min. The primer pair RpS2 
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Table I . List of specimens used for molecular analyses including GenBank accession numbers. 



Species 


Sex 


Locality and date 


Lab code 


COI 

accession 
number 


RpS2 

accession 

number 


subgenus Colias 
Fabricius, 1807 


Colias hyale (Linnaeus, 1758) 
irkutskana Stauder, 1 923 


male 


Russia, SW Transbaikalia, 
Buryatia, Selenga river district, 
Gusinoye Ozero village env., 
steppe rivulet valley, 7.6.2003 


MZHJL35 


HE775142 


HE775198 


Colias hyale (Linnaeus, 1758) 
irkutskana Stauder, 1 923 


male 


Russia, SW Transbaikalia, 
Buryatia, Selenga river district, 
Gusinoye Ozero village env., 
steppe rivulet valley, 7.6.2003 


MZHJL44 


HE775143 


HE775199 


subgenus Eriocolias 
Berger, 1986 


Colias adelaidae adelaidae 
Verhulst, 1991 


male 


China, Gansu, Xia-He, 3400 m, 
35°1TN, 102°31'E, 25.6.2004 


MZHJL61 


HE775187 


HE775243 


Colias alpherakii alpherakii 
Staudinger, 1882 


female 


Kyrgyzstan, Alai mts., 4 km 
SE Tengizbai pass, 3400 m, 
3.7.2001 


MZHJL37 


HE775169 


HE775225 


Colias alpherakii alpherakii 
Staudinger, 1882 


female 


Kyrgyzstan, Alai mts., 4 km 
SE Tengizbai pass, 3400 m, 
3.7.2001 


MZHJL51 


HE775180 


HE775236 


/""> !• 1 11 1 II T" 

Lolias berylla berylla Fawcett, 
1904 


male 


China, S Tibet, Himalaya 
Mts., Lablungla pass, 4800 m, 
18-22.7.2001 


MZHJL48 


HE775178 


HE775234 


Colias berylla berylla Fawcett, 
1904 


male 


China, Tibet, Lhodak, 4600 m, 
15.7.2002 


MZHJL55 


HE775182 


HE775238 


Colias christophi christophi 
Grum Grshimailo, 1885 


C 1 

remale 


Tadjikistan, Turkestanskyi Mts., 
Kumbel pass, 3000 m, July 2002 


MZH _JL45 




I ir ■ / cool 

HE775231 


Colias christophi helialaica 

C L 1^. i no o 

icnulte, 19oS 


male 


Kyrgyzstan, Alai Mts., W end 
of Tengizbai pass, 3700 m, 
5-6.7.2001 


MZHJL67 


HE775192 


HE775246 


Colias cocandica cocandica 
Erschoff, 1874 


male 


Kyrgyzstan, Suusamyr Mt. r., 
Alabel pass, 3200 m, 10.7.2002 


MZHJL43 


HE775174 


HE775230 


Colias cocandica hinducucica 
Verity, 1911 


male 


Tajikistan, E Pamir, Ak-Buura 
Mts., 4250 m,l4-15.7.2003 


MZHJL34 


HE775168 


HE775224 


Colias cocandica pljushtchi 
Verhulst, 2000 


male 


Kyrgyzstan, Sary Dzhaz riv. bas., 
Kaindy-Ketta mts., Tashkoro 
village, 3000 m 10.7.2003 


MZHJL19 


HE775160 


HE775216 


Colias eogene C. et R. Felder, 
[1865] elissa Grum Grshimailo, 
1890 


male 


Kyrgyzstan, W end of Tengizbai 
pass, 3700 m, 5-6.7.2001 


MZHJL1 


HE775144 


HE775200 


Colias eogene C. et R. Felder, 
[1865] elissa Grum Grshimailo, 
1890 


male 


Kyrgyzstan, W end of Tengizbai 
pass, 3700 m, 5-6.7.2001 


MZHJL40 


HE775171 


HE775227 


Colias fieldii Menetries, 1855 
chinensis Verity, 1909 


male 


China, Sichuan, Zhangia, 3000 
m, 32°47'N ) 103°36'E, 6.6.2002 


MZHJL50 


HE775179 


HE775235 


Colias fieldii Menetries, 1855 
chinensis Verity, 1909 


female 


China, Gansu, Shin-Long-Shan, 
2800 m, 35°48'N, 103°59'E, 
29.6.2004 


MZHJL60 


HE775186 


HE775242 


Colias grumi grumi Alpheraky, 
1897 


female 


China, Gansu, Altun Shan, road 
from Aksay to Danjing pass, 
2500-2800 m, 22-23.7.2002 


MZHJL54 


HE775197 
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Species 


Sex 


Locality and date 


Lab code 


COI 

accession 
number 


RpS2 

accession 

number 


Colias heos heos 
(Herbst, 1792) 


male 


Russia, SW Transbaikalia, 
Buryatia, Selenga river district, 
Gusinoye Ozero village env., 
steppe rivulet valley, 1.7.2003 


MZHJL39 


HE775170 


HE775226 


Colias heos heos (Herbst, 1792) 


male 


Russia, SW Transbaikalia, 
Buryatia, Selenga river district, 
Gusinoye Ozero village env., 
steppe rivulet valley, 1.7.2003 


MZHJL46 


HE775176 


HE775232 


Colias lada lada Grum 
Grshimailo, 1891 


male 


China, Sichuan, Maningano 
surr., 31°56'N, 99°12'E, 4500 
m, 15.6.2002 


MZHJL7 


HE775150 


HE775206 


Colias lada lada Grum 
Grshimailo, 1891 


male 


China, Sichuan, Maningano 
surr., 31°56'N, 99 < T2'E, 4500 
m, 15.6.2002 


MZHJL27 


HE775165 


HE775221 


Colias ladakensis Felder, 1865 
seitzi oollow, 1939 


male 


China, SW Tibet, Himalaya 
Mts., 100km WParyang, 4650- 
5000 m, 13.6.2004 


MZHJL4 


HE775147 


HE775203 


Lotias ladakensis welder, 186;> 
seitzi Bollow, 1939 


male 


China, SW Tibet, Himalaya 
Mts., 100km W Paryang, 4650- 
5000 m, 13.6.2004 


MZHJL57 


HE775183 


HE775239 


Colias marcopolo marcopolo 
Grum Grshimailo, 1888 


male 


Tadjikistan, E Pamir, Dunkeldyk 
Lake, 4400 m, 25.7.2003 


MZHJL30 


HE775166 


HE775222 


Colias marcopolo marcopolo 
Grum Grshimailo, 1888 


1 

male 


Tadjikistan, E Pamir, Dunkeldyk 
Lake, 4400 m, 25.7.2003 


MZH _JLii 


T_TTT"7"7C 1 / — 7 

Ht7/3l67 


OC77C ni 

HE7/3223 


Colias marcopolo marcopolo 
Grum Grshimailo, 1888 


male 


Tadjikistan, E Pamir, Dunkeldyk 
Lake, 4400 m, 25.7.2003 


MZHJL41 


HE775172 


HE775228 


Colias montium montium 
Oberthur, 1886 


male 


China, Sichuan, Maningano 
surr., 31°55'N, 99°12'E, 4000 
m, 9-18.6.2004 


MZHJL59 


HE775185 


HE775241 


Colias nebulosa Oberthur, 1894 
sungpani Bang-Haas, 1 927 


male 


China, Sichuan, Maningano 
surr., 31°56'N, 99°12'E, 4500 
m, 15.6.2002 


MZHJL9 


HE775152 


HE775208 


Colias nebulosa Oberthur, 1894 
sungpani Bang-Haas, 1927 


male 


China, Sichuan, Maningano 
surr., 31°56'N, 99 < T2'E, 4500 
m, 15.6.2002 


MZHJL24 


HE775162 


HE775218 


Colias nebulosa Oberthur, 1894 
sungpani Bang-Haas, 1 927 


male 


China, Sichuan, Maningano 
surr., 31°56'N, 99°12'E, 4500 
m, 15.6.2002 


MZHJL26 


HE775164 


HE775220 


Colias nina Fawcett, 1904 
hingstoni Kiley, 192 3 


male 


China, SW Tibet, Himalaya 
Mts., 60 km S Saga, 4600-5000 
m, 7-8.6.2004 


MZHJL53 


HE775181 


HE775237 


Lotias mna rawcett, 1904 
hingstoni Riley, 1923 


male 


China, SW Tibet, Himalaya 
Mts., Lablongla pass, 4800 m, 
5.6.2004 


MZHJL58 


HE775184 


HE775240 


Colias regia regia Grum 
Grshimailo, 1887 


1 

male 


Kyrgyzstan, Kaindy-Ketta Mt. r., 
Kumar pass, 3200 m, 12.7.2003 


A /I '/ I_T TT 0 

MZri _JLo 


HE7/31D1 


HE/732U7 


Colias regia regia Grum 
Grshimailo, 1887 


male 


Kyrgyzstan, Kaindy-Ketta Mt. r., 
Kumar pass, 3200 m, 12.7.2003 


MZHJL42 


HE775173 


HE775229 


Colias romanovi romanovi 
Grum Grshimailo, 1885 


male 


Kyrgyzstan, Alai mts., 4 km 
SE Tengizbai pass, 3400 m, 
7-8.7.2001 


MZHJL3 


HE775146 


HE775202 
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Species 


Sex 


Locality and date 


Lab code 


COI 

accession 
number 


RpS2 

accession 

number 


Lotias romanovi romanovi 
Grum Grshimailo, 1885 


male 


Kyrgyzstan, Alai mts., 4 km 
SE Tengizbai pass, 3400 m, 
7-8.7.2001 


MZHJL47 


HE775177 


HE775233 


Colias sieversi sieversi Grum 
Grshimailo, 1887 


male 


Tadjikistan, Peter I Mts., 
Ganishob, 2400 m, 17.6.2004 


MZHJL70 


HE775195 




Colias sifanica sifanica Grum 
Grshimailo, 1891 


i 

male 


China, Gansu, Xia-He, 3400 m, 
SSTl'N, 102°31'E, 25.6.2004 


A A '/ I_T TT 1 1 

MZht |L1 1 


T_TTT"7"7C 1 C / 

HE773154 


T irnr -Tin 


Colias sifanica sifanica Grum 
Grshimailo, 1891 


male 


China, Gansu, Xia-He, 3400 m, 
35°H'N, 102°31'E, 25.6.2004 


MZHJL64 


HE775189 


HE775245 


Colias staudingeri Alpheraky, 
1881 pamira Grum 
Grshimailo, 1890 


male 


Kyrgyzstan, Zaalaisky (Transalai) 
Mts., Altyn Dara river, 3000 m, 
25.7.2000 


MZHJL2 


HE775145 


HE775201 


Colias staudingeri Alpheraky, 
1881 pamira Grum 
Grshimailo, 1890 


male 


Kyrgyzstan, Zaalaisky (Transalai) 
Mts., Altyn Dara river, 3000 m, 
25.7.2000 


MZHJL13 


HE775156 


HE775212 


Colias staudingeri Alpheraky, 
1881 pamira Grum 
Grshimailo, 1890 


male 


Kyrgyzstan, Zaalaisky (Transalai) 
Mts., Altyn Dara river, 3000 m, 
25.7.2000 


MZHJL23 


HE775161 


HE775217 


Colias stoliczkana stoliczkana 
Moore, 1882 


male 


India, Jammu Kashmir, Ladakh 
Range, Markha Valley, Ganda 
Pass, 4600 m, 12.7.2001 


MZHJL15 


HE775158 


HE775214 


Colias thisoa Menetries, 
1832 aeolides Grum 
Grshimailo, 1890 


male 


Kyrgyzstan, Sary Dzhaz riv. bas., 
Kaindy-Ketta mts., Tashkoro 
village, 3000 m, 10.7.2003 


MZHJL10 


HE775153 


HE775209 


Colias thisoa Menetries, 
1 832 aeolides Grum 
Grshimailo, 1890 


female 


Kyrgyzstan, Sary Dzhaz riv. bas., 
Kaindy-Ketta mts., Tashkoro 
village, 3000 m, 10.7.2003 


MZHJL17 


HE775159 


HE775215 


Colias thisoa Menetries, 
1832 aeolides Grum 
Grshimailo, 1890 


female 


Kyrgyzstan, Sary Dzhaz riv. bas., 
Kaindy-Ketta mts., Tashkoro 
village, 3000 m, 10.7.2003 


MZHJL25 


HE775163 


HE775219 


Colias thrasibulus thrasibulus 
Fruhstorfer, 1910 


i 

male 


China, W Tibet, Mandhata Mt., 
4900 m, 15-16.7.2003 


A T '/ I_f TT 1 / 

MZri_JL14 


HE773 07 


Lime n i q 

HE775213 


Colias tibetana tibetana Riley, 
1922 


male 


China, Tibet, Himalaya Mts., 
Nyalam, 4200 m, 8.7.2003 


MZHJL6 


HE775149 


HE775205 


Colias tibetana tibetana Riley, 
1922 


male 


China, SW Tibet, Himalaya 
Mts., Nyalam, 3700^200 m, 
28-30.6.2004 


MZHJL63 


HE775188 


HE775244 


Colias wanda wanda Grum 
Grshimailo, 190/ 


male 


China, Qinghai, 20km NW 
of Zhidoi City, 4700-5000 m, 
16.7.2000 


MZHJL66 


HE775191 




Colias wanda wanda Grum 
Grshimailo, 1907 


male 


China, S. Tibet, Cona, 4500- 
4700 m, 24-25.6.2004 


MZH _JL69 


HE775194 




Colias wiskotti Staudinger, 
1882 draconis Grum 
Grshimailo, 1891 


male 


Uzbekistan, Chandalas Mts., 
Chakmksh village, 2600 m, 
27.6.2004 


MZHJL71 


HE775196 




Colias wiskotti Staudinger, 1882 
hofmannorum Eckweiler, 2000 


male 


Iran, Khorasan, 75km SE of 
Birjand, 2200 m, 18-20.5.2002 


MZHJL68 


HE775193 




Colias wiskotti Staudinger, 
1882 separata Grum 
Grshimailo, 1888 


male 


Kyrgyzstan, Alai mts., 4km 
SE Tengizbai pass, 3400 m, 
3.7.2001 


MZHJL65 


HE775190 





subgenus Eucolias 
Berger, 1986 
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Species 


Sex 


Locality and date 


Lab code 


COI 

accession 
number 


RpS2 

accession 

number 


Colias tyche tyche (de Boeber, 
1812) 


male 


Russia, East Siberia, Lake 
Baikal, Khamar-Daban Mrs., 
Slyudyanka river, taiga, 800 m, 
14.6.2003 


MZHJL5 


HE775148 


HE775204 


Colias tyche tyche (de Boeber, 
1812) 


male 


Russia, East Sayan, Buryatia, 
Mondy env., Huruma river, 
1500 m, 6.6.2002 


MZHJL12 


HE775155 


HE775211 



Table 2. List of Colias GenBank samples of the COI barcode used in this study. 



Species 


GenBank accession number 


Colias alpherakii 


FJ663407 


Colias christophi 


FJ663409 


Colias chrysotheme elena 


FJ663410 


Colias chrysotheme elena 


FJ663411 


Colias croceus 


EF457737 


Colias croceus 


FJ663412 


Colias croceus 


GU688507 


Colias croceus 


HQ004279 


Colias croceus 


HQ004282 


Colias eogene 


FJ663415 


Colias eogene 


FJ663416 


Colias erate amdensis 


EF457736 


Colias erate poliographus 


EF457735 


Colias erate poliographus 


EU583852 


Colias erate poliographus 


GU372561 


Colias fieldii 


EF584859 


Colias hyale 


FJ663418 


Colias hyale 


FJ663421 


Colias hyale 


HQ004297 


Colias hyperborea 


EF457739 


Colias marcopolo 


FJ663422 


Colias marcopolo 


FJ663423 


Colias myrmidone 


HQ004303 


Colias phicomone 


HM393178 


Colias regia 


FJ663427 


Colias tamerlana mongola 


FJ663424 


Colias tamerlana mongola 


FJ663425 


Colias tamerlana mongola 


FJ663426 


Colias thisoa thisoa 


FJ663429 


Colias tyche 


FJ663430 


Colias wiskotti chrysoptera 


FJ663431 


Colias wiskotti chrysoptera 


FJ663432 


Colias wiskotti chrysoptera 


FJ663433 


Colias wiskotti wiskotti 


FJ663435 


Colias wiskotti wiskotti 


FJ663436 
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nF (5 '-ATCWCGYGGTGGYGATAGAG-3 ') and RpS2 nR ( 5 ' -ATGRG GCTTKC- 
CRATCTTGT-3') (Wahlberg and Wheat 2008) was used to amplify a ca. 400 bp 
fragment of the nuclear RpS2 gene. The PCR were carried out following the PCR 
cycling profile described in Wahlberg and Wheat (2008): initial heating 95 °C for 7 
min, 40 cycles of 95 °C for 30 s, 50 °C for 30 s, 72 °C for 2 min, and a final exten- 
sion period of 72 °C for 10 min. Sequencing of the double-stranded PCR product was 
carried out on an ABI PRISM* 377 Automated Sequencer (Applied Biosystems) fol- 
lowing manufacturer's recommendations. All PCR primers were used for sequencing. 
Sequences were inspected and edited using Sequence Navigator* (Applied Biosystems). 



Sequence analysis 

We analysed and clustered our sequence data using parsimony and Neighbour-Joining 
(NJ) of K2P-distances. We used parsimony and NJ for our newly generated COI se- 
quence dataset, NJ for RpS2 sequences, parsimony for the concatenated COI and 
RpS2 sequences, and, finally, NJ for the combined COI sequences generated in this 
study and those in GB. All trees were rooted using Papilio glaucus (family Papilionidae) 
and Aporia crategi (Pieridae, subfamily Pierinae) as outgroup taxa. 

Parsimony analysis was performed using NONA (Goloboff 1 999) and spawn with 
the aid of Winclada (Nixon 2002), using a heuristic search algorithm with 1000 ran- 
dom addition replicates (mult* 1000), holding 10 trees per round (hold/ 10), max trees 
set to 10,000 and applying TBR branch swapping. All base positions were treated as 
equally weighted characters. Nodal support was assessed with bootstrap resampling 
(1000 replicates) using Winclada (Nixon 2002). MEGA5 (Tamura et al. 2011) was 
used for NJ clustering using 1000 bootstrap replicates. The Kimura 2-parameter mod- 
el was used for NJ clustering of the COI sequences, while the Tamura-Nei model with 
gamma distributed rates was chosen for the RpS2 sequences. 



Results 

Sequences 

We obtained a 643 bp COI barcode for 56 Colias specimens, and a 409 bp fragment 
of RpS2 was obtained for 49 specimens (Table 1). A+T content of the COI sequences 
was 69.22%, and of the RpS2 45.0%. There were 115 parsimony informative sites for 
COI and 39 for RpS2. 

Uncorrected pairwise divergences between ingroup taxa ranged between 1.09 and 
4.09% (mean 2.77%) for COI and 0.0-1.7% (mean 1.0%) for RpS2. GenBank ac- 
cession numbers are given in Table 1 . Intraspecific uncorrected distances were up to 
1.09% (in C. thisod) for COI, with specimens of most species differing by less than 4 
nucleotide changes. 
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Identification: COI vs. RpS2 

The parsimony analysis of the new COI sequences yielded four equally parsimonious 
trees (CI = 0.59, RI = 0.75) the strict consensus tree of which is presented in Figure 1. 
The NJ tree is presented in Figure 2. 

The majority of the species could be identified with COI alone, as no COI 
haplotypes were shared between species. Both parsimony and NJ trees recovered 25 
(out of 28) species as monophyletic groups (Figures 1-2). Neither Colias cocandica, 
nor C. nebulosa formed monophyletic entities, as their sequences were scattered over 
various parts of the trees. The two samples of C. tyche were not recovered as sister 
taxa, for sample MZHJL5 appeared as sister taxon of C. heos. The overall topolo- 
gies of the parsimony and NJ trees were identical, except for the placement of C. 
thrasibulus. Parsimony placed the taxon as sister to a clade of five taxa (Figure 1), 
while NJ placed it as sister to C. romanovi (Figure 2). The external morphology of 
C. thrasibulus is rather different from that of C. romanovi, while some similarities 
can be found between C. thrasibulus and C. nina, C. ladakensis, C. tibetana and C. 
cocandica (Figure 1). Only 17 of the 39 parsimony informative sites of RpS2 were 
variable among the 49 ingroup members. NJ only recovered few species as separate 
lineages due to the shallow divergences (Figure 3). The information content of this 
gene region is best interpreted as a character-based diagnostic table, as suggested by 
DeSalle et al. (2005). This gene region yielded species specific (diagnostic) haplo- 
types for 1 1 species out of 33 (Table 3). 



Analysis of the concatenated COI + RpS2 data 

The parsimony analysis of COI + RpS2 yielded nine trees of length 560 steps (CI = 
0.63, RI = 0.72), the strict consensus tree of which is shown in Figure 4. Colias cocan- 
dica, C. nebulosa and C. tyche were not monophyletic and C. thrasibulus had the same 
position as in the COI cladogram (Figure 1). 



Analysis of all the COI sequences 

The strict consensus cladogram for all the available COI data resolved the taxa in the 
same positions as in the tree of the new COI sequences only. For ten species of the 
present study sequences were also available from GB. Sequences of most species clus- 
tered together as monophyletic entities, except for C nebulosa, C cocandica, C. tyche 
and C. regia. For C. regia the GB sequence-(GB accession no FJ663427) did not cluster 
together with our sequences. The GB barcodes of C. erate and C. croceus were shared 
by these two taxa. 

Neither the Himalayan and south Tibetan adjacent mountain Colias fauna 
{berylla, ladakensis, nina, stoliczkana, thrasibulus, tibetana), nor the east Tibetan, 
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Table 3. Species haplotypes for 17 variable posirions of RpS2 for Cenrral Asian Colias species (RpS2 
data matrix positions no 14, 152, 170, 176, 189, 191, 194, 195, 218, 284, 287, 302, 341, 353, 356, 
365, 380). 



Haplotype 


positions of RpS2 


MZHJL35_hyale 


TCCCCGGGTCCATTTTC 


MZHJL44_hyale 


TCCCCGGGTCCATTTTC 


MZHJL02 staudingeri 


TCCTCGAGTTCAAATCC 


MZHJL13 staudingeri 


TCCTCGAGTTCAAATCC 


MZHJL23_staudingeri 


TCCTCGAGTTCAAATCC 


MZHJL43 cocandica cocandica 


TCCCCGAGTTCAAATCC 


MZHJL41 marcopolo 


TACCCGAGTTCAAAACC 


MZHJL30_marcopolo 


TACCCGAGTTCAAAACC 


MZHJL07_lada 


TCCCAAAAGTCGATTCC 


MZHJL27_lada 


TCCCAAAAGTCGATTCC 


MZHJL25_thisoa 


TCCCAAAAGTCGATTCC 


MZHJL10_thisoa 


TCCCAAAAGTCGATTCC 


MZHJL17_thisoa 


TCCCAAAAGTCGATTCC 


MZHJL05_tyche 


TCCCAAAAGTCGATTCC 


MZHJL12_tyche 


TCCCAAAAGTCGTTTCC 


MZHJL39_heos 


TCCCAAAAGTCGATTCC 


MZHJL46_heos 


TCCCAAAAGTCGATTCC 


MZHJL53_nina 


TCCCAAAAGTCGATTCC 


MZHJL58_nina 


CCCCCGAAGTCGATTCC 


MZH _JL1 l_sifanica 


TCCCCGAGGTCGWTTCC 


MZHJL64_sifanica 


TCTCCGAGGTCGATTCC 


MZHJL57_ladakensis 


TCCCCGAGGTCGATTCC 


MZHJL06_tibetana 


TCCTCGAGGTTATTTCC 


MZHJL09_nebulosa 


TCCTCGAGGTTATTTCC 


MZHJL26_nebulosa 


TCCTCGAGGTTATTTCC 


MZH _JL1 4_thrasibulus 


TCCTCGAGGTTATTTCC 


MZH TL01 eogene 

—J — O 


TCCTCGAGGTTATTTCT 


MZHJL04_ladakensis 


TCTCCGAGGTTATTTCC 


MZH _JL1 5_stoliczkana 


TCTCCGAGGTTGTTTCT 


MZH TL19 cocandica nliushtchi 


TCCTCGAGTTCATTTCC 


MZH_JL34 cocandica hinducucia 


TCCTCGAGTTCATTTCC 


MZHJL03_romanovi 


TCCTCGAGTTCATTTCC 


MZH JL08 regia 

—J — b 


TCCCCGAGTTCATTTCT 


MZHJL42_regia 


TCCCCGAGTTCATTTCT 


MZHJL47_romanovi 


CCCTCGAGTTCATTTCC 


MZHJL51_alpherakii 


TCCCCGAGTTCATTTCC 


MZHJL37_alpherakii 


CACCCGAGTTCATTTCC 


MZHJL67_christophi_christophi 


TCCTCGAGTTCATTTCC 


MZH JL4 5_christophi_kali 


TCCTCGAGTTCGTTTCC 


MZHJL40_eogene 


TCCTCGAGGTTGTTTCT 


MZHJL24_nebulosa 


TCCTCGAGGTCGTTTCC 
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Haplotype 


positions of RpS2 


MZH JL5 9_montium 


CCCTCGAGGTTGTTTCC 


MZHJL61_adelaidae 


TCCTCGAGGTCGTTTCC 


MZHJL60_fieldii 


TCCTCGAGGTTATTTCC 


MZHJL50_fieldii 


TCCTCGAGGTTATTTCT 


MZHJL33_marcopolo 


TCCCCGAGGTCATTACT 


MZHJL63_tibetana 


TCCTCGAGGTTATWTCC 


MZHJL48_berylla 


TCCCCGAGGTCGAATCC 


MZHJL55_berylla 


TCCCCGAGGTCGAATCC 



- PAPILIO GLAUCUS 



- APORI A CRATAEG1 




JLUthrasibulus 
JL53 nina 
JL58 nina 

JL61 

JL6 tibetana 
JL63 tibetana 
[— o — JL19 cocandica pljushtchi 
\ ~ L-oo — JL34 cocandica hinducucia 

"^-oooo-I - JU l8dak8nsis 
1 — JL57 ladakensis 
JL3 romanovi 
JL47 romanovi 

JL24 nebulosa 

JL43 cocandica cocandica 
JL10 thisoa 

JL17thlsoa 
JL25 thisoa 

JLiSstoliczkana 

JL66 wanda 
JL69 wan da 

J L1 3 staudingerl 

JL2 staudingerl 
JL23 staudingerl 
JL70 sieversi 

JLB regia 
JL42 regia 
JL1 eogene 
JL40 eogene 
JL41 marcopolo 

JL30 marcopolo 
JL33 marcopolo 
JL45 chrlstophi kali 
JL67 chrlstophi chrlstophi 

JL27 lada 

-oo— j |-<xx>*- JL9 nebulosa 

| JL11 slfanlea 

JLG4 sifanica 
JL12tyche 

JL5 tyche 
JL39 heos 
JL46 heos 
JL37 alpherakii 
JL51 alpherakii 

JL68 wiskotti hofmannorum 
JL65 wiskotti 
JL71 wiskotti draconis 



Figure I . Strict consensus cladogram of Colias COI sequences obtained in this study. 



Qinghai, Gansu and Sichuan species aggregates {adelaidae, grumi, lada, montium, 
nebulosa, sifanica, wanda) were resolved as species clusters similar to the Tian Shan, 
Pamir and Hindukush species. 

Several COI haplotypes were noted for a few species, even among specimens ob- 
tained from the same locality (e.g. C. staudingeri and C. thisoa), Taxa not resolved 
as monophyletic clusters were the species C. cocandica and C. nebulosa. All the in- 
cluded subspecies of C, cocandica (G c. cocandica, C. c. pljutshtshi and C c. hinducucia) 
showed distinct COI sequences, with cocandica cocandica as most different. 
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JL65_wiskotti 

JL68_wiskotti_hofmannorum 
JL71__wiskotti_draconis 
JL37_alpherakii 
94 L JL51„alpherakii 
r JL12_tyche 

Hi JL5_tyche 

4 jj JL39_heos 
9? JL46_heos 

r JL45_christophi_kali 
94 JL67_christophi_christophi 
99 | JL7_lada 

\- JL27_lada 
— JL9_nebulosa 
;gh JL11_sifanica 
85^ JL64_sifanica 

JL43__cocandica_cocandica 
JL10_thisoa 

JL17_thisoa 
99' JL25_thisoa 

JL24_nebulosa 

JL15_stoliczkana 
J — JL66_wanda 

JL69_wanda 
■ JL70_sieversi 
llj" JLIeogene 
JL40_eogene 
98 r JL8_regia 
JL42_regia 
r JL41_marcopolo 
^"[| JL30_marcopolo 
77 JL33_marcopolo 
i JL13_staudingeri 
JL2_staudingeri 
75 ' JL23_staudingeri 

i JL14_thrasibulus 

3] v~ JL3_romanovi 

99 1 Jl_47_romanovi 
100 1 JL6_tibetana 
' JL63_tibetana 

' JL61_adelaidae 

84 [ JL53_nina 
JL58_nina 
99 i JL4_ladakensis 
JL57Jadakensis 

) r JL19_cocandica_pljushtchi 

99 JL34_cocandica_hinducucia 
i JL59_montium 
~v\ i JL50 fieldii 



100' JL60_fieldii 
JL26_nebulosa 

Jl_54_grumi 
JL48_berylla 
100' JL55_berylla 
JL35_hyale 



100' JL44_hyale 



PAPILIO_GLAUCUS 



APORIA_CRATAEGI 



Figure 2. Neighbour-Joining tree using the K2P-parameter model for the COI sequences obtained in 
this study. 
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r! 



JL39_heos 
JL46_heos 
JL12_tyche 
JL5_tyche 
JL17_thisoa 
JL25_thisoa 
JL10_thisoa 
JL27_lada 
JL7_lada 
JL53_nina 
JL58_iiina 
JL11_sifanica 
JL57_ladakensis 
JL48_berylla 
JL55_berylla 
JL64_sifanica 
JL4_ladakensis 
^6 1 JL15_stoliczkana 
JL24_nebulosa 
JL61_adelaidae 
- JL59_montium 

JL34_cocandica_hinducucia 
JL3_romanovi 
JL1 9_cocandica_pljushtchi 
JL23_staudingeri 
JL13_staudingeri 
JL2_staudingeri 
~ JL67_christophi_christophi 
" JL47_romanovi 
■ J L45 ch ri stophi kal i 
50J JL8_regia 
JL42_regia 
JL33_marcopolo 
83j JL41_marcopolo 
JL30_marcopolo 
JL37_a!pherakii 
95 1 JL35_hyale 
L JL44_hyale 
1 S| — JL43_cocandica_cocandica 

JL51_alpherakii 
JL6_tibetana 
JL9_nebulosa 
JL26_jiebulosa 
JL14_thrasibulus 
JL60_fieldii 
JL63_tibetana 
■ JL40_eogene 
33 li JLI eogene 
W JL50_fieldii 

APORIA_CRATAEGI 

PAPILIO_GLAUCUS 



0.02 

Figure 3. Neighbour-Joining tree using the Tamura-Nei model with gamma distributed rates for the 
RpS2 sequences. 
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Figure 4. Strict consensus cladogram of the concatenated data set of COI + RpS2. 



Discussion 

Barcoding 

Lukhtanov et al. (2009) tested the utility of COI barcodes for Central Asian butterflies 
by sampling specimens from a considerable geographical range. They observed that 
this substantially increased intraspeciflc variation reducing the interspecific divergences 
("barcoding gap"), but that this did not hamper species identification. The present 
study shows that most Colias taxa form monophyletic entities that can be identified 
with COI data alone. The RpS2 gene region showed identical sequences in cocandica 
pljutshtshi and cocandica hinducucia (Table 3, Figure 3), differing by only three nucleo- 
tides from cocandica cocandica. Based on the molecular data the recognition of these 
subspecies is not or weakly supported. 

The fact that the three C. nebulosa samples were scattered over different parts of the 
COI tree might be the result of a laboratory contamination due to carry over between 
samples. The C. nebulosa samples were collected on the same day and in the same 
place. C. nebulosa is morphologically distinct from other Colias species, excluding pos- 
sible misidentification. The RpS2 data, however, could point to two morphologically 
cryptic species in sympatry (samples MZHJT24 vs. MZHJT9 and MZHJT26), so 
that the different COI barcodes might represent numts, despite no apparent 'signs' (no 
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indels). This discrepancy between morphology and DNA sequence data emphasises 
the necessity to use multiple samples to detect this sort of challenging issues. 

Even though C. cocandica and C. nebulosa did not form monophyletic groups 
our results show that COI barcodes are useful for (1) identifying Palaearctic and 
Central Asian Colias, (2) pointing to a possible cryptic species, and (3) highlight- 
ing the necessity to further investigate the question on the subspecific rank of C. 
cocandica cocandica. 

The utility of RpS2 as a species barcode for Colias spp. is clearly more limited, since 
e.g. C. heos, C. lada, C. nina, C thisoa of the subgenus Eriocolias and C. tyche (subgenus 
Eucolias) have identical sequences (Table 3, Figure 3). Still, RpS2 yielded species spe- 
cific (diagnostic) haplotypes for 1 1 species of the subgenus Eriocolias and for C. hyale 
(subgenus Colias s.str.). 



Congruence with traditional classification: analysis of concatenated COI + RpS2 

The strict consensus tree was more resolved than either of the trees resulting from sepa- 
rate analyses of the gene regions (Figure 4). 

Although the concatenated data did not resolve the phylogenetic relationships 
among all Colias species, some observations can be made. The majority of the species 
confined to the adjacent Tian Shan, Pamir and Hindukush mountain ranges form 
a well supported clade. This includes C. eogene, C regia, C. romanovi, C. marcopolo, 
C staudingeri, C. christophi, C. alpherakii and C. wiskotti. Yet, C sieversi, which also 
occurs in these mountain ranges (Peter I and Khozratishoh mountains), was not in- 
cluded in this clade. C. sieversi is morphologically most similar to C. alpherakii, thus 
showing another case of disagreement between morphological and DNA sequence 
data. C. thisoa, too, lives in the aforementioned mountain ranges, but it has a wider 
distribution, stretching from Turkey to the Altai Mountains. A third taxon, C. c. co- 
candica, is considered closely related to C. tamerlana (e.g. Verhulst 2000), a species 
occurring in southern Siberia and Mongolia. Thus, the origin of C. thisoa and C. c. 
cocandica may differ from that of the species confined to the Tian Shan, Pamir and 
Hindukush mountain range. One sample of C. cocandica (MZHJL43) was placed 
within this "mountain" clade, while the other two samples appeared as sister taxa to 
the Himalayan species C. ladakensis. As with C. sieversi, our DNA data disagree with 
the morphological characters, but it should be noted that this clade is not well sup- 
ported. Conversely, two morphologically similar Himalayan species, viz. C. nina and 
C. ladakensis, were assigned to different clades. In the COI + RpS2 tree they were 
placed in different, more encompassing species clusters (Figure 4), in the COI NJ tree 
they were joined with C. c. pljutshtshi and C. c. hinducucia (Figure 2), while the COI 
cladogram resolved these taxa together with C. adelaidae, C. tibetana, C. c. pljutshtshi 
and C. c. hinducucia (Figure 1). 
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The analyses did not support the monophyly of the subgenera Eucolias and Erio- 
colias sensu Berger (1986). The Eucolias species C. tyche was not resolved as a separate 
monophyletic lineage, but was resolved into Eriocolias. This is congruent with the results 
of Pollock et al. (1998) and Brunton (1998). Only the the subgenus Colias, here repre- 
sented by C. hyale, is supported as a distinct lineage, placed as sister to all other Colias sp. 



Barcodes of Palaearctic Colias spp. 

The parsimony (Figure 5) and NJ analyses (Figure 6) of the larger matrix of Palaearctic 
COI barcodes (total COI) recovered the same species clusters, but some of the species 
show different placements (e.g. C. thisoa, C christopht). This is not surprising as all in- 
ternal nodes are very shallow. The samples of C. tyche and C. hyperborea show very low 
sequence difference, morphologically these taxa are different, and they largely share the 
same distribution area. An example of species that share the same distribution and that 
exhibit clear morphological similarities, and which as such were resolved as sister spe- 
cies in both analyses, includes C. wiskotti and C. alpherakii. Identification of Palaearctic 
Colias based on COI barcodes is in most cases possible, since shared haplotypes were 
recorded only for C. erate and C. croceus. 

Intraspecific variation is notable between some of the recognized subspecies, both 
among our own samples and those downloaded from GB. The intraspecific variation can 
partly be explained by morphologically clearly distinct subspecies, such as those of C. 
wiskotti, or by specimens from widely different localities, such the different specimens of 
C. hyale (sample FJ663418 from Russia, FJ663421 from Kazakhstan, HQ004297 from 
Romania and MZHJL35 and MZHJL44 from SW Transbaikalia). However, notable 
intraspecific variation also occurs within populations, such as C. thisoa aeolides with all 
samples originating from the same locality and date, but the limited sampling prevents 
conclusions on the reasons for this. It is apparent that the understanding of intraspecific 
variability of the COI barcode for Colias is presently very limited. 

The combined COI data of our sequences and sequences downloaded from GB 
include species belonging to one additional subgenus, Neocolias, represented by C. 
myrmidone and C. erate. Only the subgenus Colias, represented by C. hyale, is well 
supported as distinct lineage. Yet, one specimen of C. hyale (FJ6634 1 9) clustered to- 
gether with C. erate {Neocolias) and C. croceus {Eriocolias). The other subgenera were 
not resolved as clades according to present classification, in agreement with our results 
for the combined analysis. 

Our findings generally support COI as a species specific barcode for Colias, but we 
also highlight the necessity of including multiple individuals of species in molecular 
barcoding studies. Problematic 'cases' of widely divergent barcodes or conflicting mor- 
phological and molecular 'signals' are found in most if not all barcoding studies, and 
this study makes no exception. 
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Figure 5. Strict consensus cladogram of COI sequences for Palaearctic Colias taxa. 
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Figure 6. Neighbour-Joining tree using the K2P-model of COI sequences for Palaearctc Colias taxa. 
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